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I.  ANALYSIS 

In  this  section,  the  stiffness  and  compliance  transformations  are 
computed  for  the  following  stress  sign  convention: 
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Tne  computation  is  carried  out  in  a  stepwise  fashion  to  illustrate  the  points 
of  similarity  of  the  two  properties.  The  stress  rotational  transformation  is 
given  by  the  following  matrix: 


m2-n2 


where 


m  =  cos  e 
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The  preceding  matrix  also  applies  to  strain  if  y/2  is  used  for  the  shear  term 
in  the  strain  tensor: 
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For  the  present  purpose,  it  is  more  convenient  to  use  a  slight  modification 
of  T  to  define  a  transformation  matrix  for  engineering  strain: 
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It  is  noted  that  the  engineering  strain  transformation,  which  uses  t  instead 
of  >'2,  is  the  negative  transpose  of  the  conventional  stress  transformation 
matrix.  Tne  basic  mono-piy  properties  in  the  natural  directions  are  as 
follows : 


E^  ~  Longitudinal  Young's  modulus 
E-j  -  Transverse  Young's  modulus 

-  Shear  modulus 

A? 

-  Major  Poisson  ratio  defined  as  -  - —  for  load  along  L 
>  -  u  j  E  l 
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The  compliance  and  stiffness  forms  of  Hooke  s  Law  for  stress  in  the  ply 
natural  directions  are  as  follows: 
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In  order  to  use  Hooke's  law  in  the  p:y  natural  direction,  the  general  stresses 
and  strains  must  be  transformed  to  the  ply  natural  directions.  First,  we 
review  the  compliance  form: 
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Starting  with  Hooke's  law 


transforming  strain  and  stress 


and  solving  for  strain 
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Equation  (1)  is  the  compliance  form  of  Hooke's  law  for  the  anisotropic 
plate.  A  similar  analysis  of  the  stiffness  form  of  Hooke's  law  gives 
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Transforming  the  stress  and  strain  gives 
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and  solving  for  the  stress  yields  the  following  matrix: 
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The  preceding  equation  is  the  stiffness  form  of  Hooke's  law  for  the 
anisotropic  plate.  The  detailed  forms  of  the  compliance  and  stiffness 
matrices  are  shown  below  for  reference: 
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Stiffness  Matrix 
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Many  authors  [for  example,  Tsai  (Ref.  1)]  have  shown  compact  forms  for  the 
rotationally  transformed  stiffness  matrix,  using  multiple-angle  trigonometric 
identities  and  convenient  coeffients  called  U.  The  same  approach  is  taken 
here  to  simplify  the  rotationally  transformed  compliance  matrix.  In  this 
case,  the  compliance  coefficients  are  called  W.  Although  the  individual 
components  have  some  formal  similarity,  there  are  notable  detailed 
differences,  as  may  be  seer,  in  the  matrix  elements  and  coefficients  for 
stiffness  and  compliance  given  below: 


Stiffness 

Matrix  Elements 

Q( 1 1 )  =  U( 1 )+Ul2)cos(2e)+U(3)cos(46) 
Q( 22 )  =  U( 1 )+U(2)cos(26)+U(3)cos(4e) 
Qi 12)  =  U(4)  -U(3)cos(4e) 

Q(66)  =  U( 5 )  -U(3)cos(4e) 

Q(16)  =  sin(26)  +U( 3  )s in ( 46 ) 

Q(  26  )  =  sin(26)  -U(3)sin(4o) 


The  forms  of  the  preceding  matrix  elements  are  slightly  different  for  elements  with 
index  66 ,  16 ,  and  26 . 
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Coefficients 


Coeff icients 


U(1)  =  1 /8[ 3Q( 11 )+3Q(22)+2Q( 12 )+4Q(66 ) ]  W(1)  =  1/8[ 3S( 11 )+3S(22)+2S( 12)+S(66)  ] 


U(2)  =  1/8[4Q( 11 )-4Q(22) ] 

U( 3)  =  1/8[Q( 1 1 )+Q(22)-2Q( 12)-4Q(66) ] 
U( 4 )  =  1/8[Q( 1 1 )+Q(22)+6Q( 12) -4Q( 66 ) ] 
U(5)  =  1/8[Q( 1 1 )+Q(22)-2Q( 12 )+4Q(66) ] 
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W(3)  =  1/8[S( 11 )+S(22 )-2S( 12) -S(66) ] 
W(4 )  =  1 /8[ S( 11 )+S(22)+6S( 12)-S(66) ] 
W(5 )  =  1/8[S( 11 )+S(22)-2S( 12)+S(66) ] 


Also  W(5)  = 


W(1)  -  W ( 4 ) 
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The  forms  of  the  preceding  coefficients  are  slightly  different  for 
coefficients  with  index  3,  4,  and  5. 
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II.  THE  TWO  QUASI-ISOTROPIC  ASYMPTOTES 


A.  STIFFNESS- ISOTROPIC  ASYMPTOTE 

The  parallel  model,  or  stiffness-isotropic  asymptote,  corresponds  to 
uniform  random  orientation  of  plies  in  a  contiguous  laminate: 

2  71 

Q(iso)  =  J  Q( i, j,e)de 
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This  well 

(Ref.  3)] 

known  result  [see,  for  example  Ashton 

is  given  by  tne  isotropic  Hooke's  law 

et  al . 
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(Ref.  2)  and 

shown  below: 
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Tne  above  isotropic  stiffness  matrix  a. so  results  from  certain  lamination 
patterns  suer,  dS  0  sou  and  0  90  r45.  In  terms  of  the  mono-ply  engineering 
parameters  ir.  tne  natural  direction,  tnese  stiffness  matrix  elements  are  given 
DV 
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The  stiffness-isotropic  engineering  parameters  are  given  by 


v  =  U  ( 4 )  /  U  (  1) 


G  =  U(5) 


E  r  U( 1)  ( 1  -  v*)  (5) 


B.  COMPLIANCE-ISOTROPIC  ASYMPTOTE 

In  this  case  we  integrate  compliance: 

S(iso)  =  J  S(i,j,6)d8 
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The  resulting 
similar  but  no 
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In  terms  of  the  mono-ply  engineering  parameters  in  the  natural  directions, 
these  compliar.ce-matrix  elements  are  given  by 
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W  ( 4  ) 

W(5 ) 


(7; 


and  the  compliance-isotropic  engineering  parameters  are  given  by 


Direct  formulas  for'  the  pseudo-isotropic  moduli,  in  terms  of  the  ply 
natural  direction  values,  are  given  below. 
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The  aDove  equations  were  usea  in  TaPie  1  to  compute  the  two  sets  of 
quasi-isotropic  properties  for  typical  fiber  composite  materials  used  in 
aerospace  structures.  These  properties  are  analogous  to  the  Voigt  (parallel) 
and  heuss  (series  i  moaels. 
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III.  APPLICATIONS 


Analysis  of  stress  and  strain  in  rectangular  orthotropic  materials  has 
uncovered  a  number  of  geometrically  axisymmetric  problems,  in  which  the 
resultant  stress  is  axisymmetric,  despite  the  orthotropic  material  properties. 
The  corresponding  strains  in  the  cylindrical  coordinate  system  are  not 
axisymmetric,  and  of  course,  in  the  design  critical  material  natural 
directions,  neither  stresses  nor  strains  are  axisymmetric. 

The  compliance  isotropic  material  analysis  given  in  the  previous  sections 
turns  out  to  be  embedded  in  this  class  of  problems.  The  strain  energy  in  a 
ring  shaped  element  for  conventional  2D  axisymmetric  FE  analysis  is  given  by 
(see,  for  example,  Refs.  A  and  5) 
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For  an  axisymmetric  geometry  dVL  -  rd6dA1 
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where 
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Equation  (11)  can  be  recast  in  terms  of  the  stresses  by  substituting  the 


compliance  form  of  Hooke's  law  e.  =  S. ,o.. 
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This  formulation  is  a  slight  variation  of  the  line  followed  by  Paraoen 
(Refs.  4  and  5),  who  attempted  to  rationalize  a  "weighted  average"  elasticity 
or  stiffness  matrix  for  use  in  axisymmetric  finite-element  analysis  of 
rectangular  orthotropic  material.  Pardoen  cast  the  problem  in  terms  of  the 
stiffness  matrix  and  the  strains,  neither  of  which  is  axisymmetric. 

The  case  of  axisymmetric  stress  allows  Eq.  (11)  to  be  rewritten  as 


rA  _  2r 

U .  =  — T—  o.  [J  [S .  .  ]d6 I  0  . 
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The  term  in  Drackets  is  the  quasi -isotropic  compliance  matrix  derived  in 
Section  B,  and,  for  a  given  orthotropic  material,  produces  the  quasi -isotropic 
compliance  properties  as  specifically  defined  in  the  matrix  of  Eq.  (6). 

Further  stress  analysis  with  this  strain  energy  formulation  is  isotropic  and 
axisymmetr ic .  We  have,  in  effect,  defined  a  fictitious  or  equivalent 
isotropic  material  with  compliance-isotropic  properties.  Further  analysis  can 
proceed  with  the  use  of  the  corresponding  Ec,  Gc,  and  vc  of  Eq .  (8)  in  closed 
form  isotropic  solutions  where  the  stress  is  axisymmetric.  Examples  are  given 
below  for  disks  or  annul  i  of  rectangular  orthotropic  material  subjected  to 
inertial,  me char. i Ca. .  arc  thermal  loads.  These  were  used  as  test  cases  by 
Pardoen  (Refs.  U  and  5). 

Example  1:  Rotating  Disk  with  Rectangular  Orthotropy 

Stress  analysis  of  the  rotating  disk  of  radius  F  witn  rectangular 
orthotropy  (Refs.  6  and  7)  gives: 

Radial  stress-orthotropic  analysis 

op  =  (1  -  6)  ( 1  -  (1U) 

R  2  Rc 

Circumferential  (Hoop)  stress-orthotropic  analysis 
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Stresses  in  an  isotropic  disk  (e.g.,  Ref.  8)  are 


Radial  stress-isotropic 


=  e?s!  ( 3-1-2)  (,  -  4; 


Hoop  stress- isotropi c 
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1  ♦  3v  r_- 
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The  form  of  Eqs.  (19)-(’5)  is  identical  to  Eqs.  (1?)-(18)  if  1  -  6  = 


3  -  v 


With  the  definition  of  £  given  by  Eq .  (16),  it  turns  out  that  this  is  true, 
and  the,  "effective"  Poissor.  Ratio,  v,  is  the  compliance  isotropic  value 
vc  given  in  Eq  .  ( 9c)  . 

Thus  the  orthotropic  factor  s  is  a  disguise  for  a  pseudo- isot ropi c 
material  with  compliance  isotropic  properties. 

Example  2:  Pressurized  Ortnotropic  Disk 

The  rectangular  orthotropic  pressurized  annulus  has  an  axisymmetric 
stress  state,  according  to  Lekhnitskii  (Ref.  9),  if  the  orthotropic  material 
properties  satisfy  the  following  relation: 


=2-  -  2v  -  =2  s  i 
G  xy  E 

xy  y 


This  may  be  rewritten  as 


,  *  1  =  Ui-  -  2. 

s  Gxy  ^ 


Substituting  this  constraint  in  Eqs.  (9a),  (9b),  and  (9c)  for  each 
engineering  parameter  Ec ,  Gc,  and  reveals  that  this  axisymmetric  solution 
corresponds  with  an  equivalent  compliance  isotropic  material  having 


G  =  G  , 
c  xy 


Note  that  this  compliance  isotropic  property  set  defined  by  the  above 
constraint  is  explicitly  independent  of  the  transverse  modulus. 

The  radial  and  circumferential  stresses  are  given  by  the  well  known 
Lame  equation  for  thick  wall  annulus  (e.g.,  Ref.  8). 


Example  3:  Thermal  Load  on  Orthotropic  Disk 

An  axisymmetric  stress  distribution  occurs  in  a  rectangular-orthotropic 

p  p 

disk  subject  to  axisymmetric  temperature  of  the  form  T(r)  =  TQ  -  £T(r  /a  ), 

and  the  solution  is  determined  for  the  case  a  =  a  by  use  of  the  compliance 

x  y 

isotropic  modulus  in  the  equations  for  radial  and  circumferential  stress: 
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In  each  of  the  above  examples,  the  stress  is  axisymnetric  and  may  be 
computed  by  standard  isotropic  formulas  provided  the  appropriate  compliance- 
isotropic  properties  from  Eq .  (9)  are  used.  The  strains  are  not 
axisynmetric .  The  nonuniform  strain  distribution  is  obtained  by  the 
orthotropic  Hooke's  law  of  Eq.  (1). 

The  compliance  integration  and  corresponding  isotropic  asymptote  are 
found  in  analysis  of  wrinkle  defects  in  composite  materials  (Ref.  9). 


IV.  CONCLUSIONS 
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Tr.e  second  isotropic  asymptote  of  planar  fiber  composite  materials, 
derived  from  compliance  summation,  has  been  developed  here,  with  rotational 
matrix  parameters  cast  in  multiple  angle  form.  The  analytical  form  closely 
resembles,  but  is  not  identical  to,  the  stiffness  isotropic  matrix  elements 
and  coefficients. 

Formulas  for  the  isotropic  engineering  parameters  are  given  in  terms  of 
the  orthotropic  mono-ply  engineering  parameters. 

Tne  compi lance- isotropi c  properties  turn  up  m  stress  analysis  of 
orthotropic  materia,  problems  where  axisymmetric  stress  occurs.  Several  of 
tr.ese  problems  are  given  here.  The  approach  here  snows  how  this  class  of 
problems  can  ce  modeled  and  analyzed  as  isotropic  materials,  often  with 
simple,  closed  form  equations.  This  approach  might  be  used  as  an 
approximation  in  problems  wnere  ax isymme tr ic  stress  can  be  assumed. 
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